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ABSTRACT: The purpose of this short paper is to revisit the infinite Lie group theoretical frame- 
work of hydrodynamics developped by V. Arnold in 1966. First of all, we extend this approach from 
the Lie pseudogroup of volume preserving transformations to an arbitrary Lie pseudogroup. Then 
we prove that, contrary to what could be believed from the work of Arnold which is of a purely 
analytical nature, the same results can be obtained from a purely formal point of view. Finally, we 
provide the analogue for both the so-called "body" and "space" dynamical equations. We conclude 
by showing that even this new approach can be superseded by dynamics on Lie groupoids, along 
ideas pioneered by the brothers E. and F. Cosserat or H. Weyl, on the condition to change the 
underlying philosophy. 

KEY WORDS: Dynamics, Hydrodynamics, Lie group, Lie pseudogroup. Lie groupoid, gauge 
theory, Spencer operator, differential sequence, adjoint operator, duality theory. 

INTRODUCTION: 

In a celebrated paper published in 1966 [1], V. Arnold applies the differential geometry of 
infinite dimensional Lie groups to the hydrodynamics of perfect uncompressible fluids. His tech- 
nique amounts to use for the Lie pseudogroup of volume preserving diffeomorphisms of a bounded 
Riemannian domain, the analytical analogue of the formal methods already developped by H. 
Poincare [13], N.G. Chetaev [5,6], E. and F. Cosserat [7,8,9,11], G. Birkhoff [3] and himself [1,2] 
for the dynamics of a rigid body and, more generally, for the dynamics that can be achieved on 
any Lie group. It must be noticed that the last reference does not quote the previous ones and 
we advise the reader to compare [3, p 205 to 216, in particular the formula on p 215/216] with 
[2,Appendice 2, in particular the formula of Th 2.1 on p 326]. 

The basic tool, which is crucially used, is the so-called " orthogonal decomposition theorem" of H. 
Weyl roughly saying that any vector field on the previous domain can be decomposed uniquely 
into the sum of a divergence-free vector field on this domain, tangent to the boundary, and the 
gradient of a univalent function. 

Since the time we read this paper in the seventies, we have always been convinced that this an- 
alytical approach could be replaced by a purey formal approach following, for Lie pseudogroups, 
the one already existing for Lie groups, the latter one being the founding stone of gauge theory 
[21]. It is only now that we have been able to succeed while revisiting once more Arnold's paper 
in the light of new results recently obtained for the partial differential (PD) optimal control theory 
and the corresponding multivariable variational problem with differential constraints [17]. 
We first present a new interpretation of the dynamics on Lie groups in the form of a differential 
sequence called gauge sequence and exhibit the corresponding linearized sequence. 
We then extend the preceding results to Lie pseudogroups and obtain the corresponding differential 
sequence both with the corresponding variational calculus with constraints. All these results are 
new and the particular case of a one-parameter (time) gauging of the Lie pseudogroup of volume 
preserving transformations provides the equations of Arnold plus the corresponding " space" equa- 
tions that he did not obtain. The corresponding second section has been deliberately written in a 
self-contained way though it relies on tricky calculations involving the implicit function theory. 
Finally, revisiting the proof of one key theorem within the previous specific case of hydrodynamics. 
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we provide doubts about the usefulness of this approach and hints for using Lie groupoids and Lie 
algebroids instead of Lie pseudogroups and Lie algebras of vector fields (For a nice introduction to 
Lie groupoids, see [12, appendix]). 

In the case of Lie group actions, studying only the linear framework for simplicity while introduc- 
ing the Spencer operator, we exhibit the link existing between the linear gauge sequence and the 
corresponding Spencer sequence for Lie algebroids [14,15,16,18]. But this is just the generalized 
Cosserat theory that we already published [14]. 

As global actions on manifolds may not exist, all the results in this paper are local ones though 
global notations are used for simplicity in order to avoid using explicit open sets and charts. 

DYNAMICS ON LIE GROUPS: 

This section, which is a summary of results already obtained in [14,15], is provided for fixing 
the notations and the techniques leading to the gauge sequence and the corresponding variational 
calculus, both with the respective linearized versions. 

Let X be a manifold of dimension n with local coordinates x = {x^,...,x"') and latin indices 
i,j — l,...,n. A point on X will play the part of n parameters. We denote by T = T{X) the 
tangent bundle to X and by T* = T*{X) the cotangent bundle to X while AT* is the bundle 
of r- forms on X. Let now G be a Lie group of dimension p with identity e, local coordinates 
a = (a^, flP) and greek indices p, a, r. We denote hy Q ^ Te{G) the corresponding Lie algebra 
with vectors denoted by greek letters A, /i, z/. As usual, we shall identify a map a : X G called 
gauging of G over X, with its graph X ^ X x G which is a section of a trivial principal bundle, 
and, similarly, use the same notation for a bundle and its sheaf of (local) sections as the back- 
ground will always tell the right choice. In particular, when differential operators are involved, the 
sectional point of view must automatically be used. Such a convention allows to greatly simplify 
the notations at the expense of a slight abuse of language. 

If we have a map a : X — > G : a; — > a{x), we obtain the tangent map T{a) : T{X) — > Ta{G) : dx — > 
da = ^dx and we can pull back the image to Te{G) = Q hy acting with the inverse a~^ of a, either 
on the left to get A = a~^da e T* (g) ^ or on the right to get B = daa~^ & T* ^Q. Differentiating 
a~^a = e, we get ar^da = —ad{a~^) and thus daa~^ = —ad{a~^), that is one obtains B from A by 
changing a to a~^ and changing the sign too. Also we obtain symbolically A = a~^Ba = Ad{a)B 
by introducing the adjoint map fj, = Ad{a)X obtained by carrying X & Q from e to a on the left 
and coming back to € Q from a to e on the right. For more technical details on the adjoint map, 
we refer the reader to [15,p 180]. 

Another way to present the previous construction is to introduce the Maurer-Cartan left invari- 
ant 1-form uj = {uo'^{a)da'^) on G with value in [Q) and pull it back on X by T(a) in order 
to get A = AJdx'^ = ui'^{a{x))dia'^ {x)dx^ . Finally, using the Maurer-Cartan structure equa- 
tions for ui on G, namely dw'^ — c^^uj'' A co'^ =0, and pulling them back on X similarly, we 
get a well defined operator T* (g) Q —^ A^T* (g) Q : A ^ dA - [A, A] = F with local coordinates 
F[j = djA'j — djAJ — Cp^A^Aj where the c are (care to the sign) the structure constants on G 
with Lie algebra bracket ([A, ^j)"^ = c^^XPii". We may collect these results in the following theorem: 

THEOREM 1: In the previous framework, there exists the so-called gauge sequence: 

XxG ^ T* ®g A^T* (g) g 

a a~^da = A 

A ^ dA-[A,A] = F 

REMARK 1: The previous results can be extended to connections on principal bundles and their 
curvature [15] but it is important to notice that, in both cases, the group is not acting on the base 
space. 

REMARK 2: When n = 1, no differential sequence is existing and this is the situation con- 
sidered by Arnold for the only "time" parameter t. In the case of a rigid body moving in M.^, 
changing slightly the notations with Xq the initial position and x the final position at time t, the 
movement of the rigid body is a; = a{t)xo + b{t) and the projection of the speed in the body is 
a~^x = a~^dxo + a~^b with standard notation for time derivative. This result brings out at once 
the form A appearing for the Lie group of rigid motions where G = {a,b) with parameters a for 
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rotations and b for translations. Nevertheless, when n = 2, with parameters the curvilinear abcissa 
s and the time t, the above framework is the one adapted to the Kirchoff theory of thin elastic 
beams along the work of E. and F. Cosserat for the group of rigid motions [9]. Also, when n = 3, 
the two previous operators are exactly described by the brothers Cosserat in the nice reference [11] 
which has never been quoted elsewhere. An important but tricky question raised by mechanicians 
was thus to understand why the compatibility conditions of Cosserat theory were first order PD 
eqations while they were known to be second order PD equations in classical elasticity theory [11]. 

REMARK 3: The second operator has been introduced by E. Cartan for introducing the so- 
called curvature (G is the rotation group) and torsion (G is the translation group), but with no 
reference to the first operator [4]. It must be noticed that the "abstract" group G is not acting 
on X. Therefore, if one wants to relate the above framework to electromagnetism (EM), the only 
possibility is to consider X as space-time and to call A the 4-EM potential, F the EM field (made 
by E and B combined together), with the necessary condition to have dim{Q) = 1. This was the 
birth of gauge theory with G = U{1), the unit circle in the complex plane and G the parallel to 
the complex axis at he point (1,0). 

It just remains to introduce the previous results into a variational framework. For this, taking 
into account Remark 3, we need to consider a lagrangian on T*®Q, that is an action W = ^ w{A)dx 
where dx = dx^ A ... A dx" and to vary it. We obtain successively; 

5A = 5a~^da + a~^6da 

= —{a~^da){a~^da) + a~^d{aa~^Sa) 

= d{a~^6a) + {a~^da){a~^6a) — {a~^da){a~^da) 

Setting a-^Sa = A G = A°T* O Q, we thus obtain [13,14,15]: 

dA = d\- [A, A] 
Finally, setting dw/dA = A= (Al) G A^-^T* Q, we get: 

6W = j A6Adx = j A{dX - [A, \])dx 

and therefore, after integration by part, the Euler-Lagrange (EL) equations (with no right mem- 
bers) [1(13), 15]: 

We notice that such a linear operatorfor A has non constant coefBcients linearly depending on A. 
However, setting 5aa~^ = /x, we also have: 

SA = —a~^Saa~^da + a~^d((Saa~^)a) 

= —a~^6aa~^da + a~^d{Saa~^)a + a~^Saa~^da 
= Ad{a)dfj. 

Therefore, introducing by duality B such that Bfj, = AX, we get the equivalent form [1 (14), 15]: 

diBl = 

which is a divergence-like operator. 

We let the reader check by himself, as an exercise, the following other formulas for right invariant 
objects: 

6B = diJ,+ [B, fj] = Ad{a-^)dX 
where A is used in place of fj, and vice- versa, while is used in place of a and signs are changed. 

Similarly, setting A ""^^ = ada~^ = —B and caring to the sign, we obtain therefore 

dB + [B, B] = 0. This formula will be found again later on in a different framework. 
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At the end of this section, we provide the "infinitesimal" Hncar version of the previous results 
for a "close" to e. Then there is no difference between A and B or between A and fi, and we get 
the linear gauge sequence: 

A°r* ®g ^ A^T* ®g ^ a^t* ® g 

which is just the tensor product by ^ of a part of the Poincare sequence for the exterior derivative. 

DYNAMICS ON LIE PSEUDOGROUPS: 

Before using specific notations for this particular section in a way coherent with the first section, 
let us provide an elementary introduction to the local theory of Lie pseudogroups with notations 
that will also be used in the third section where they will also be coherent with the notations of 
the first section too. 

DEFINITION 1: A Lie group of transformations of a manifold X is a lie group G with an action 
of G on X better defined by its graph X x G ^ X x X : {x,a) ^ {x,y = ax = f{x, a)) with the 
properties that a{bx) = {ab)x and ex = x,Va; € X,ya,b € G. 

It is sometimes useful to distinguish the source x from the target y by introducing a copy Y oi X 
with local coordinates y = (y^, y"). Such groups of transformations have first been studied by S. 
Lie in 1880. Among basic examples when n = 1 we may quote the affine group y = ax + b and the 
projective group y = [ax + h) / {cx + d) of transformations of the real line. When n = 3 we may quote 
the group of rigid motions y = ax + b where now a is an orthogonal 3x3 matrix and 6 is a vector. 
Such a group is known to preserve the euclidean metric u = {lo^^ = u)^^) and thus the quadratic 
form ds'^ = {dx^)'^ + (dx^)^ + (dx^)"^ = uj'^^dx^dx^ . When n = 4 we may quote the conformal group 
of space-time with 15 parameters (4 translations, 6 rotations, 1 dilatation, 4 elations) preserving 
the Minkowski metric lo or the quadratic form ds^ = (dx^)"^ + (dx^)^ + (dx^)"^ — c^{dt)'^ up to a 
function factor, where now c is the speed of light and t the time. Among the subgroups, we may 
consider the Weyl group with 11 parameters preserving uj up to a constant factor and the Poincare 
group with 10 parameters preserving lo. We recall that the three Lorentz transformations must be 
considered as space-time rotations, the three other rotations being pure space rotations. 

Only ten years later, in 1890, S. Lie discovered that the Lie groups of transformations were 
only examples of a wider class of groups of transformations, first called infinite groups but now 
called Lie pseudogroups. 

DEFINITION 2: A Lie pscudogroup F of transformations of a manifold X is a group of transfor- 
mations y = f{x) solutions of a (in general nonlinear) system of partial differential (PD) equations, 
also called system of finite Lie equations. 

Roughly, this definition means that, ifUq{X, Y) with local coordinates {x,yq) — (x, y'',yf,yfj, ...) 
satisfying det{yf) 7^ is the q-jet bundle of invertible transformations (just replace derivatives by 
symbols !), there is a system TZq C Uq{X,Y) defined by equations $'^(a;,t/g) = such that, if we 
have two solutions that can be composed, the composition is again a solution. However, such a point 
of view cannot be tested in actual practice. Instead, the idea is to use sections of TZq, that is maps 
fq : (x) ^ (x, /'^(x), //=(x), /^-(x), ...) satisfying det{ff{x)) ^ such that $-(x, /,(x) = 0, Vx G X. 
Introducing the specific section jq{f) : (x) {x,f''{x),dif''{x),dijf''{x),...) and the composi- 
tion jq{g o f) = jq{g) o jq{f) for maps that can be composed, an equivalent definition that can 
be tested is that, whenever fq and gq are sections of TZq that can be composed, then gq o fq is 
also a section. A similar comment can be done for defining the inverse f~^ and is left to the reader. 

Setting now y = x + t(_{x) + ... and passing to the limit for t — > 0, that is to say linearizing TZq 
around the q-jct of the identity y — x, we get a linear system Rq C Jq{T) for vector fields with solu- 
tions © C r satisfying [O, Q] C O. It can be proved, for the same testing type reasons, that such a 
system may be endowed with a Lie algebra bracket on sections ^q : (x) — > (x,^'^(x),^*^(x),^*^-(x), ...) 
that we shall quickly define (see [14,15] for more details and compare to [12]). Such a bracket 
on sections transforms Rq into a Lie algehroid and we have [Rq,Rq] C Rq. Let us first de- 
fine by bilinearity {jg+i(^), jg+i(7y)} = ig([^, ??]), V^, € T. Introducing the Spencer operator 
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D : Rg+i -^T* (g)Rg : ^g+i jii^g) - ^g+i with local components {di^'' - - ^fj, ...), we 

obtain the following general formula at order q: 

[£.q^Vq] = {£.q+l,Vq+l} + KODllg+l - i{T])D^g+l ,\/£,g, r]g £ Rg 

where iQ is the interior multiplication of a 1-form by a vector, and we let the reader check that 
such a definition no longer depends on the "lifts" (,g+i,i]g+i over £,q,r]q. 

EXAMPLE 1: (Affine transformations) n=l,q = 2,X = R^ 

With evident notations, the system TZ2 is defined by the single linear OD equation y^x = and 
the sections are defined by fxx{x) = 0. Similarly, the sections of i?2 are defined by txxi^) = 0. Ac- 
cordingly, the components of [^2, il^] at order zero, one and two are defined by the totally unusual 
successive formulas: 

v] = S.dxV - vdx^ 

{[^2,V2])xx = ixTlxx - rjxixx + ^dxTJxx - VdxS.xx 

It follows that i^xx = 0, r]xx = ([^2, »72])a;x = and thus [i?2, -R2] C i?2- 
EXAMPLE 2: (Projective transformations) n = 1, g = 3, X = 

The system 7^3 is defined by the single nonlinear OD equation [yxxx/Vx) ~ ^{yxx/Vx)^ — and 
the sections of i?3 are defined by S,xxx{x) — 0. The formulas for the bracket of Lie algebroid 
[i?3, -R3] C -R3 can be derived similarly but involve many more terms. 

EXAMPLE 3: (Volume preserving transformations) n arbitrary, q = 1,X = M" 

The system TZi is defined by the single nonhnear PD equation d{y^ , y^)/d{x^, x") = det{yf) = 

1 and the sections of Ri are defined by the single relation Q = 0. Accordingly, we obtain: 

([6, r?l])? = ^Iv'r - nl^r + CdrV! ' dr^^ 

When summing on k and i, the first two terms disappear (as in Example 1 !) and we get therefore 
[Ri,Ri] C R\. We invite the reader to compare this result with the usual way on solutions where 
one defines O as the kernel of the Lie derivative C{£)io of the (volume) n-form w = dx^ A ... A dx^ 
with respect to ^ and then uses the well known formula [C{^) , C{r])]u) = £([^, r/])a;, V^, e T in 
order to obtain [0, 6] C 0. 

REMARK 4: In view of the results of the first section, it is tempting to consider a Lie pseu- 
dogroup of transformations as an infinite Lie group of transformations. Such a point of view, 
frequently adopted in the past by russian mathematicians (V. Arnold, L.V. Ovsiannikov,...) is at 
the opposite of point of view adopted by the american school (D.C. Spencer and coworkers). The 
origin of this confusion lies in the fact that the transformations of a few Lie pseudogroups just 
depend on a certain number of arbitrary functions. Expanding these functions into Taylor series 
close to the value they have at the identity transformation allows to consider the coefficients of 
the series as parameters on an infinite Lie group. The best way is by far to obtain results from the 
finite number of defining finite or infinitesimal Lie equations, using their sections without using 
their solutions. The following tricky example will illustrate this remark. 

EXAMPLE 4: Let us consider the Lie pseudogroup of transformations of only depending on 
an arbitrary function of a single variable with nonzero derivative: 

r = {2/ = fix ),y =a;/-^-j— } 
It is easy to check that F is the set of solutions of the (involutive) nonlinear system: 

^1 y2-o,yyi-x ^^^^^^^^^ -i 



5 



with corresponding system of infinitesimal Lie equations: 



Ri c Ji(T) ^2 = 0, x^ei +e =0^^1+^1=0 

We may thus set f{x) = ao + a\x + ... with a\^Q and / must be the identity f{x) = x whenever 
y = X ,y = X . 

Nevertheless, changing slightly the last system to the following one: 



does not allow now to have any generic explicit solution. We finally notice that, in the first case 

one can check the Lie pseiidogroup property by composing explicit solutions while in the second 
case the corresponding Lie pseudogroup is the one preserving the geometric object uj = (a, /?) where 
a = x^dx"^ — x^dx^ is a 1-form and /? = dx^ A dx^ is a 2-form satisfying the integrability condition 
da = 2/3. 

Prom now on in this section, let X be the manifold of parameters as in the first section, let Y 

be another manifold of dimension m with coordinates y = (y^, y™) and indices k, I while Z is a 
copy of Y with coordinates z = {z^, z™). Let now F = Ty Q aut{Y) be a Lie pseudogroup of 
transformations of Y {not of X) of the form z = f{y) with inverse y = g{z). With a slight abuse 
of language, a map / : X ^ F : x — > fx{y) = f{v,x) will be identified with a section of X x F. 
We denote by Rq = RqiY) C Jq{T(Y)) the corresponding system of infinitesimal Lie equations. 
The operator T> = Vy : T{Y) — > Jq{T (Y)) / Rq{Y) is a Lie operator, that is an operator such that 
V^i = 0, 2?^2 = => C2] = 0, where the bracket is the standard bracket of vector fields on Y. 
Again with a slight abuse of language, if 6 = Qy C T{Y) is the sheaf of germs of solutions of V, 
we have [6, 6] C 6. 

The main problem that we have now to solve is to construct an analogue of the gauge se- 
quence with left and right points of view, both with a corresponding variational calculus. For this, 

considering fx+dx °fx^ — id and passing to the limit, we obtain a vector valued 1-form (CARE) 
f € T* (g) Qz CT* (Si V{X X Z) called generalized speed and defined by the formula: 

v{z,x) = ^{g{z,x),x) 

with no one of the indices for simplicity. Such a definition is just similar to the definition of the 

eulerian speed v(x, t) in continuum mechanics for a movement x = /(xq, t) of a point/particle with 
position xq at time to ; a; at time t and x + v{x, t)dt at time t -\- dt. However, the Lie pseudogroup 
background is not so evident and the reader may now, as a motivation, get in mind the composition 

z ^ y — > z + v{z, x)dx. 

Similarly, considering /^T^ o fx+dx — id and passing to the limit, we may obtain the pull-back 
u G T* O Or C T* (g) V{X X Y) defined by the formula: 

df df 

-g^{y,x) = v{f{y,x),x) = —{y,x)u{y,x) 

again with no one of the indices for simplicity as only the tangent mapping is involved. 
Conversely, if y = g{z,x),y = g{z,x) € F, then, from the implicit function theorem, we get 
y = h{y, a;) e F in general. However, from the identities: 

/(.9((^, x), x) = z = f{g{z, x), x) , f{h{y, x),x) = f{y, x) , g{z, x) = h{g{z, x), x) 
we deduce the relations: 

dy dx ' dy dx 



and thus: 



dq dhdq dh ^ , , s ^ 

v{z,x) = viz,x)^- = --^- = 0^y = h{y)Gr 
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a result showing the right invariance of the generahzed speed. 

Exactly as in the previous section, we notice that the passage from v to u just amounts to change 
f to g while changing the sign (CARE) as we have indeed the identity: 

y = 9{f{y,x),x),Vx e x 

and, differentiating with respect to x, we get: 

^{z,x)v{z,x) + ^{z,x) = 

and thus: 

dg 



or equivalently: 



u{g{z,x),x) = - — {z,x),\/x G X 



dg 

u{y,x) = - — {f{y,x),x),yx e X 



but u,v € T* <Si O do not provide the same sections. 

The next important step will be to provide the compatibility conditions. For this, introducing: 

= %T(9iz,x),x) 

we successively obtain from the chain rule for derivatives and the above comment: 

dvf ^ dg' I 0^/'° 

9x3 dx'dy'- dxo dx'dxo 

dz^ dy'- dx^ dx^dx^ 

dz" ~ dx'dxi 

Exchanging i and j, then substracting, we get the formula: 

diVj - djVi + [vi,Vj] = 

that we may rewrite in the condensed form dv + [v,v] = as v can be considered as an horizontal 

as 

dx 



form on X X Z. Also, as in the first section, with u in place oi A, v in place of B and ^ in place 
of = —A, while caring about the sign, we have similarly du — [u, u] = 0. 
Collecting these results and taking into account that I> is a Lie operator with dtV = 0, Vi = 1, n, 
we obtain therefore: 

THEOREM 2: There exists a differential sequence: 

Axr — >T* (^e — > A^T* ® e 

The hard step in this paper will be to construct the corresponding variational calculus with 
constraint. For this, 6 being the usual symbol for variation, we may introduce the vertical vector 
field T] = r]{z,x) G 6 such that ri{f{y,x),x) = Sf{y,x) and we get (compare to [7,14,20]): 

dv dri dri dn . dri , . 

+ r]— = +V-!- = ^ 5v = + [v,rf\ 

oz ox oz dx ox 

Exactly as above and for a later use, by analogy with continuum mechanics, we may introduce the 
Jacobian determinant: 

dif] df 
A{z,x) = -^{g{z,x),x) = det{-^){g{z,x),x) 
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Using the well known formulas for jacobian determinants: 



d{f + df)d{f) _ d{z + vdx) _ dv 
d{y) ' d{y) d{z) dz 

and passing to the limit, we obtain therefore ^ = and thus: 

dA + f]— = A— ^ ^SA = A— ^ - 77— 

az az oz oz 

In order to overcome the minus sign in this formula and simplify it, we may therefore introduce 
the analogue p = 1/A of the mass per unit of volume and obtain the variation [5,12,18]: 

^ dp dr] - dipv'') 

We notice that such a result could not be obtained in [1] where A = l hy assumption, though this 
is coherent with the divergence-free condition for the speed in this case. 

As in the first section, we may also introduce the pull back of r] as the vertical vector ^ = 
^(2/, x) G © defined by the formula: 

df da 
r]{f{y,x),x) = —{y,x)^{y,x) ^ —{z,x)r]{z,x)=^{g{z,x),x) 

exactly as we got u from v but now with 6 in place of d. We obain the following tricky theorem, 
coherent with the results of the first section but absent from [1,2]: 

THEOREM 3: = f^H- 

Proof. Using the well known result saying that the bracket of vector fields commutes with the 
action of any diffeomorphism while working ONLY with {z, x) and CARING ABOUT THE SIGN, 
we successively get from the above formulas: 

~ af + laf ^' af ^] 

_ 9g_dri _ idg_ ^1 
dz dx '■dx^'^i 

dx\dz'll dx\dz)'l \.dx''^i 
~ dx^dydx dz^dx''! dx dy ^ dz\ dx ' 

= ^{g{z,x),x) 



as we have indeed: 



thus ending the proof. 



Q.E.D. 



By analogy with continuum mechanics, among general action integrals of the form W = 
J w{^, ^)dyAdx, we shall only consider the action integrals of the form W — J p{z, x)w{v{z, x))dzA 
dx and vary them. As in [7(V),14], the best procedure in actual practice is to introduce a "small" 
parameter e and consider the family f{y, x; e) g{z, x; e) with variation 5 = d/de also leading to 
the previous results when {z, x; e) are considered as independent variables. Of course, we have to 
add the constraint z = fx{y) = f{y,x) € F and this will be the delicate point to overcome. 
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Introducing V = ^ = (V^) or its dual U we get: 

Using X in place of e for p, we get the " conservation of mass" identity: 

Finally, integrating by part over X x Z exactly as for Eulerian coordinates (a;, t) in classical hy- 
drodynamics [7,14,20], we obtain the EL equations over X x Z, up to sign: 

or equivalently over X xY: 

-^e = 0, G e. 

It is just the solution of this problem that we have already solved for creating PD optimal control 
theory. We consider two situations: 

1) If P can be parametrized by an operator then any rj G can be written in the form 
V-i9 = T] for an arbitrary 9 G V{X x Z) in the previous variation. Then, denoting by ad the 
formal adjoint of an operator and integrating by part with respect to z, we find ad{T>-i)£ = as 
final EL equations. 

2) In the general situation, if {$^} is a fundamental generating set of differential invariants 

of r at order q and {ti)'^(y)} their value at the identity, the Lie form of the finite Lie equations 
defining F at order g is — ix>'^(y) = in the sense that we have a defining groupoid of order q. 
It follows that the constrained variational calculus for an action density w amounts to vary the 
unconstrained action density w + A^($'^ — Lu'^{y)) where {A^} arc Lagrange m,ultipliers. As the 
vertical bundle of the groupoid is isomorphic to the corresponding infinitesimal Lie equations over 
the target, changing A to A by duality, the variation becomes 5w + XDrj. Integrating by part as 
usual while caring about the signs, we obtain the EL equations £ — ad{T>)X = 0. Finally, we just 
need to eliminate the Lagrange multipliers by introducing compatibilty conditions for ad{T>) of the 
form ad{T>-\) and obtain the same EL equations as before but with a slightly different approach. 

EXAMPLE 5: When n = l,m = 3,1^ = with Euclidean metric, time t, initial position 
xo, final position x and transformations of the form x = f{xo,t), the above v{x,t) is just the 
ordinary speed u in Eulcr variables and u{xo,t) its differential pull back in Lagrange variables. If 
F = aut{Y) and the action density is the standard kinetic energy w = ■^piP where p is the mass 
per unit volume, the EL equations become, up to sign [5(V),12]: 

^^''dt^^^di^ ^'^■'^^'^^ = ^^dt + ^ ^) A iT+ V(-i?2)) = 0. 
Now, if F is the Lie pseudogroup of volume preserving transformations defined byA = l<^4'/Ci=l, 
then 6 is the Lie algebra of divergence- free vector fields and we get: 

7.77 = 0, V77 e 9 <^4> V.^ = 

In that specific case, T> = div can be parametrized by the curl operator, that is ff = ^ A for an 
arbitrary vector field 6. As the curl operator is self-adjoint up to sign, introducing the so-called 
"vortex" vector uj = 5V A v (see below for the numerical factor), we get at once the EL equations 
of Arnold, namely: 

1^ ^ div r-, ^1 

in a purely formal way. 

We point out that the numerical factor 1/2 comes from the fact that, in the case of a rigid body 
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motion x = a{t)xQ + b(t), then v = x = cixq + b = da~^x + (6 — da^^b). An easy computation in 
local coordinates then shows that da^^ is a skew-symmetric 3x3 matrix amounting to the vector 
product by the vortex vector to as previously defined. 

Otherwise, considering the action density w-|-A(A — 1) and varying it, we get Sw + X\7 .ff. Integrat- 
ing by part, we get 7 -|- VA = 0=^'VA7 = 0as before but the Lagrange multiplier A = AA = A 
now plays the part of the pressure, as well known in fluid dynamics. 

EXAMPLE 6: With m = 1, setting z' = dz/dy and so on, Le Lie form of the group of projective 
transformations of the real line is $=^^7 where $ is the well known Schwarzian 

derivative. Using the prolongation formulas 5z = rj^Sz' = rj'z' = (Sz)',... and so on, it is not so 
easy to check that (5$ = {z')'^r]"' and the needed isomorphism is thus A = {z')'^X. 

EXAMPLE 7: From everybody diving experience or looking at corks on the surface of the sea, a 
natural model for swell is well described by water particles moving along circles with radius smaller 
and smaller with depth. A parametrized form in the vertical (a;, y) plane can be: 

X = a + R{b)cos{Lot — ^(a)) , y = b + R{b)sin{ujt — (p{a)) 

As we have the well known formula S^'^'^K = ^Jf'll / ^'f,°'^^'' , the movement is that of an uncom- 

o(XQ,yQ) d(a,b) ' d(a,b) ' 

pressible fluid iff §^ = = 1 + RR'f' + {W + R')sin does not depend on time, that 

is to say ^'{a) = —R'/R{b) = k and thus <^(a) = ka + c,R{b) = Roe~'^^. In that case, it is not 
at all evident to see directly that the movement is NOT stationary in the fixed frame. However, 
for proving this result, one just needs to notice that the trajectories are bounded circles while the 
stream lines are infinite cycloidal curves (consider the surface of the sea !). Finally, the movement 
becomes stationary in the moving frame x = x — jt, y = y associated with a boat "surfing" on the 
crest of a given wave. Indeed, changing the parameters (a, 6) to (a = a — ^t, b = b), we notice that 
the new parametrization (a, b) (x, y) does not involve time t anymore. 

DYNAMICS ON LIE GROUPOIDS: 

Let us revisit the variation formula provided by Theorem 3, changing the notations when n = 1, 
one parameter t and movement x = f{xo,t) with inverse xq = g{x,t) as usual in continum me- 
chanics where the speed in Eulerian coordinates (x, t) is defined by : 

d f 

v{.f{xo,t),t) = -^{xo,t). 

Using space-time coordinates {x^ ^) with x"* = at and reference speed of light c = 1, we 

may divide {dx^ , dx^ , dx"^ , dx'^) by dt and obtain the extended speed (w^, i;^, w^, 1) = {v, 1) in Special 
Relativity with variation (6v,0). Introducing the "time" variable to corresponding to xq and the 
function to = g'^{x, t), we may finally consider the transformation x = /(xq) with inverse xq = g{x), 
but now on space-time, and try to obtain 6v = [v, 77] on space-time, on the condition to get rid 
of to- For this we refer to [8,12] as it is out of the scope of this paper. In particular, the strange 
form dx A dt is no longer a wedge product of forms but a volume form on space-time as we have 
no longer TWO separate manifolds but only ONE, a result leading to relativistic mechanics. 
Having in mind the first section where G was an abstract group, we have now a Lie pseudogroup 
of transformations of a manifold, a particularly simple example being produced by a Lie group 
of transformations defined by the graph XxG— >XxXofan action. In this case, it is known 
[14,15] that the (nonlinear) Spencer sequence is isomorphic to the gauge sequence and the Spencer 
(nonlinear) sequence is therefore the ONLY candidate for a generalization to cases where no Lie 
group action is involved, like the Lie pseudogroup of volume preserving, symplectic or contact 
transformations. 

Wc have already explained in many books [14,15] that such an aproach, where sections of Lie 
groupoids/jct bundles generalize the gauging approach, is JUST the way to understand the Cosserat 
theory of continuum mechanics and Weyl theory of electromagnetism. The only choice remaining 
on space-time is to select a convenient "group candidate" and we refer to the previous books for 
discovering why it MUST be the conformal group of space-time endowed with the Minkowski met- 
ric. 
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Finally, exactly like in the first section, a linearized version exists and we refer the reader to 
[16,17,19] for a simple presentation of this framework with applications to field-matter coupling. 



Coming back to mathematics for a few lines, there is another differential sequence to be found 
in the literature and that we did not speak about, namely the Janet sequence, though it is for sure 
the best known differential sequence. For short, if E, F, Fq, Fi, ... denote vector bundles over X, 
we use the same letters for the corresponding sets (sheaves to be exact) of sections and such an 
interpretation must be used whenever operators are involved. Starting from a vector bundle E (for 
example T) and a linear differential operator T> : E ^ F : ^ ^ rj, if we want to solve the linear 
system with second member = rj even locally, one needs " compatibility conditions" in the form 
Virj = 0. Denoting now F by Fq, we may therefore look for an operator T>i : Fq ^ F^ : rj ^ C, and 
so on. Under assumptions on I? which are out of the scope of this paper (involutivity !), the french 
mathematician M. Janet has proved in 1920 that such a chain of operators ends after n steps and 
we obtain the linear Janet sequence, namely [15]: 

^ Q ^ E ^ Fo^ Fi^ ...^ Fn ^ 

It follows that we only have at our disposal for any application where group theory seems to be 
involved, three linear differential sequences, namely the Janet sequence, the Spencer sequence and 
the gauge sequence. As these sequences are made by quite different operators, the use of one 
excludes the use of the others. 

In order to escape from this dilemna, at the end of this paper and for the sake of clarifying the 
key idea of the brothers Cosserat by using these new mathematical tools, we shall explain, in a way 
as elementary as possible while using only the linear framework, why THE JANET SEQUENCE 
AND THE GAUGE SEQUENCE CANNOT BE USED IN CONTINUUM MECHANICS. By this 
way we hope to convince the reader about the need to use another differential sequence, namely the 
SPENCER SEQUENCE, though striking it could be. Also we shall use very illuminating examples 
in order to illustrate our comments. 

First of all we exhibit the isomoprphism existing between the linear gauge sequence and the linear 
Spencer sequence. For this, if now G acts on X with a basis — {Cr^fc} of infinitesimal generators, 
we may introduce the bundle Jq{T) of g-jets of T over X, that is the vector bundle over X with 
sections transforming like the derivatives of vector fields up to order q, and the map: 

A^T* ^g^Jg{T): A^(x) ^ X^x)d^^^{x) 

It is known [14, p. 308] that this map becomes injective for q large enough and we may call Rq its 
image for such a, q. It follows from its definition that Rq ~ Rq+i is a system of infinitesimal Lie 
equations of finite typo and wo get for the Spencer operator [12,15,17]: 

D : Rq+i ^T*(S>Rq: ^,+1 ^ {Mf, - C^iJ = d,V {x)d^C';{x) 
Accordingly, the linear gauge sequence is isomorphic to the linear Spencer sequence: 

^ e ^ A°T* ^Rq^ A^T* ®Rq^ A^T* Rq 

the three isomorphisms being induced by the (local) isomorphism X x Q ^ Rq just described 
above. It is essential to notice that, though the linear Spencer sequence and the isomorphisms 
crucially depend on the action, by a kind of " miracl^^ the linear gauge sequence no longer depends 
on the action. 

This result proves that the linear Spencer sequence generalizes the linear gauge sequence, with the 
major gain that it can be used even for Lie pseudogroups of transformations that are not coming 
from Lie groups of transformations, as we shall see eamples in the sequel. 

REMARK 5: When n = 3 and we deal with the Lie group of rigid motions, the corresponding EL 

equations are nothing else but the formal adjoint of the (first) Spencer operator. Surprisingly, this 
is EXACTLY the result found by the brothers Cosserat, namely the so-called stress and couple- 
stress equations for Cosserat media [9, p 137,14,16,18,19]. 
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REMARK 6: The above result, in perfect agreement with the piczzoelcctric or photoelastic cou- 
phng of elasticity and electromagnetism, CONTRADICTS gauge theory where the lagrangians are 
functions on A^T* g and NOT on T* (S) ^ as in the previous remark [14,16]. 

Let us consider a (finite) volume jydV in limited by a (closed) surface S = j gdS and let 
us introduce the outside unit normal (pseudo) vector n = [rij) on S. Let us now suppose that the 

surface element dS is acted on by the outside with a force dF = adS and a couple dC = j7,dS, 
where both a and jl linearly depend on n through the stress tensor density a = (cr*-') and the 
couple-stress tensor density fx = (fj,^'^^ = — /i'''^'). It must be noticed that, using the standard 
Cauchy tctrahcdral device, there is no reason "a -priori" to suppose that the stress tensor is sym- 
metric. We also suppose that the volume element dV is acted on by (see later on for the sign) a 
force —fdV and a momentum —rhdV with / = (p) and m = (m*^ = — m-'*). 
Our purpose is now to study the equilibrium of the corresponding torsor fields with respect to an 
arbitrary cartesian frame Ox^x^x^. 

The equilibrium of forces is satisfied if we have the relation: 

/ adS- [ fdV = / a'^mdS - [ fdV = 

J S J V J S J V 

Using Stokes formula, this is equivalent to the well known stress equations: 

dicj'i = P 

This result shows that the surface density of forces & is equivalent, from the point of view of force 
equilibrium, to a volume density of forces / and this interpretation explains the sign adopted. 
Finally, the equilibrium of forces being satisfied, it is known that the equilibrium of momenta 
is also satisfied if it is satisfied with respect to an arbitrarily chosen cartesian frame. Hence, in- 
troducing the vector r = [x^^x^, x^), the equilibrium of momenta is satisfied if we have the relation: 

/ {il+fAa)dS- [ {rh + fAf)dV = 
J s J V 

Projecting onto the axis Ox^, we obtain: 

/ {ii^'^^+x^a^^-x'^a^^)nrdS- [ {m^^ + x^ f - x"^ f)dV = 

J S J V 

Using again Stokes formula and the previous stress equations, we obtain the couple-stress equations: 

drf/''^ + a'^ - a^' = m'^ 

This result shows that the surface density of forces a and couples jl is equivalent, from the point 
of view of torsor equilibrium, to a volume density of forces / and to a volume density of momenta 
m, provided the preceding stress and couple-stress equations are satisfied, and this interpretation 
explains the sign adopted. 

The combination of the stress AND couple-stress equations have first been exhibited by E. and F. 
Cosserat in 1909 [8,9,pl37] WITHOUT ANY STATIC EQUILIBRIUM EXPERIMENTAL BACK- 
GROUND and we now explain the key argument leading to the same equations just from group 
theoretical arguments. Of course, most of the engineering continua such as steel, concrete, glass, 
wate,... have the specific "constitutive laws" = 0,m = and we obtain therefore ct'-' = cr^% that 
is the stress tensor is symmetric, a situation not always encountered in liquid crystals. 

First of all, for the reader not familiar with the Spencer operator, we exhibit a similar result 
in a quite simpler 1-dimensional situation that will allow to recapitulate all the previous results.. 

EXAMPLE 8: Let us consider the Lie group of afiine transformations of the real line defined 

by the group action y = a^x + . The corresponding 2-dimcnsional Lie group G has coordinates 
a = {a} ^a}) and the group composition law is ah = {a^ ,a^){b^ ^b"^) — {a^b^,a^b^ +a'^) with inverse 
law = (1/a^, — a^/a^). Accordingly, we have a~^da = {{1 / a^)da^ , {1 / a^)da^) and obtain at 
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once the Maurer-Cartan forms uj^ — {1 / a^)da^ , oj'^ — {l/a^ jda'^ with the two Maurcr-Cartan equa- 
tions duj^ — 0, duj'^ + uj^ A Lo'^ — 0. Finally, a corresponding basis of infinitesimal generators of the 
action may be obtained with = x-^,^,^ = and we have in a coherent way [$^,^^] = — C^, 
that is the only non zero structure constant is = — 1- It follows from the first section that the 
resulting EL-equations are either: 

Si^l + AlAi^ = , diJ^ - AlA^ = 

or simply: 

diBi = , diB'2 = 

if we use a pidl-back by the adjoint action. In both cases these equations could not have anything 
to do with the stress and couple-stress equations previously exhibited. 

Let us now deal with the Spencer sequence instead of the gauge sequence in this framework. 
First of all, we may consider the above Lie group of transformations as a Lie pseudogroup defined 
by the second order system of finite Lie equations y^x = 0. The corresponding system R2 C J2{T) 
of infinitesimal Lie equations is ^xx = and the isomorphisms between the gauge sequence and 
the Spencer sequence is induced by the maps: 

{X\x),X^x)) ^ i^x) = xX\x) + X^{x),Ux) = XHx),^xx{x) = 0) 
The only two non-zero components of the Spencer operator become: 

dx^{x) - L{x) = xdxX^{x) + dxX^{x),dx(,x{x) - = dxix{x) = dxX^{x) 
Equating to zero these two components amounts to have: 

dxX^ =0 , dxX'' = 

Accordingly, gauging A just amounts to choose an arbitrary section of i?2- 

The final touch, that could not be in the mind of any reader even on this very simple example, is to 
work out the formal adjoint of the Spencer operator. For this, multiplying the first component by 
a test function cr(x), the second by a test function fi{x), then summing and integrating by parts, 
we get the EL-equations (up to sign) as the kernel of the following operator with second members 
(/,m): 

dx(^ = f , dxH + cr = m 

The comparison with the previous mechanical results needs no comment. 

Taking into account this example, we now study the foundation of elasticity theory and we 
restrict the study to 2-dimensional (infinitesimal) elasticity for simplicity as the general situation 
has already been treated elswhere and we just want to explain why the only founding problem of 
elasticity is the choice of an underlying Lie pseudogroup and an adapted differential sequence. 

1) The gauge sequence cannot he used: 

Looking at the book [9] written by E. and F. Cosserat, it seems at first sight that they just con- 
struct the first operator of the nonlinear gauge sequence for one parameter [9, p 7], two parameters 
[9, p 66], three parameters [9, p 123] and finally four space-time parameters [9, p 189]. This is 
NOT TRUE indeed because, according to the comment done in the first section or Example, the 
EL-equations are either a divergence like operator or a linear operator with coefficients depending 
on A, a situation not met in the couple-stress equations which is a linear operator with constant 
coefficients, not of divergence type;. In fact, a carefull study of the book proves that somewhere 
the action of the group on the space is used, but this is well hidden among many very technical 
formulas (Compare [9] p 136 with [14] p 295). 

2) The Janet sequence cannot be used: 

This result is even more striking because ALL texbooks of elasticity use it along the same scheme 
that we now describe. Indeed, after gauging the translation by defining the " displacement field^ 
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^ = (^^(x), ^^(x)) of the body, from the initial point x — {x^,x'^) to the point y = x + £,{x), 
one introduces the (small) "deformation tensor'' e = l/2£(^)a; as one half the Lie derivative with 
respect to ^ of the euclidean metric w, namely, in our case, the three components (care): 

e = (en = 01^,^12 = £21 = l/2(5i^2 ^ g^^i^^ ^ g^^2^ 

From the mathematical point of view, one uses to consider the Lie operator I?^ = £{^)u -.T ^ S2T* 
(symmetric tensors), sometimes called Killing operator, through the formula: 

One may check at once the only generating "compatibility condition" "Die = 0, namely: 

3iie22 + 922611 - 25i2ei2 = 
which is nothing else than the Riemann tensor of a metric, linearized at u. 

However, the main experimental reason for introducing the first operator of this type of Janet 
sequence is the fact that the deformation is made from the displacement and first derivatives but 
must be invariant under any rigid motion. In the general case it must therefore have {n + n^) — 
(n + n{n — l)/2) = n{n + l)/2 components, that is 3 when n = 2, and this is the reason why 
introducing the deformation tensor e. For most finite element computations, the action density 
(local free energy) to is a (in general quadratic) function of e and people use to define the stress 
by the formula cr'-' = dw/deij which is not correct because w only depends on en, £12, £22 when 
n = 2 as the deformation tensor is symmetric by construction. Finally, textbooks escape from this 
trouble by deciding that the stress should be symmetric and this is a vicious circle because we have 
proved it was not an assumption but an experimental result depending on specific constitutive 
laws. Accordingly, when n = 2, we should have (j'^-'cij = a^^eu + (2(T^^)ei2 + o-^^e22- Hence, even 
if we find the correct stress equations with this convenient duality keeping the factor " 2" , we have 
no way to get the stress AND couple-stress equations TOGETHER. 

3) Only the Spencer sequence can be used: 

Let us construct the formal adjoint of the Spencer operator by multiplying all the (2 x 2) + 2 = 6 
linearly independent nonzero components by corresponding test functions. For simplifying the 
summation, we shall raise and lower the indices by means of the (constant) euclidean metric, set- 
ting in particular — oJirSJ" and = uJir£,j- Comparing to Example 8, the only nonzero first jets 
coming from the 2x2 skewsymmetric infinitesimal rotation matrix of first jets are now ^12 = —^2,1 
while the second order jets are zero because isometrics are linear transformations. We obtain in 
the present situation: 

+ <T'^{dl^2 - 6,2) + (T^Hd2il - 6,1) + ^''52?2 + M"''9,a,2 

Integrating by parts and changing the sign, we just need to look at the coefficients of ^1 , ^2 and 
^1,2, namely: 

a — ^ Sia" + 92(721 = /I 

6,2 a^/*'-'i2 + C7I2 - C721 = mi2 

in order to get the adjoint operator ad{D) : a"~1T* (g) Rl ^ ^nj^* ^ ^* . ^^-j _^ ^-j relating 
for the first time the torsor framework to the Lie coalgebroid Rl . These equations are exactly the 
three stress and couple-stress equations of 2-dimensional elasticity. In the n-dimensional case, a 
similar calculation, left to the reader as an exercise of indices, should produce exactly the n(n+ 1)/2 
stress and couple-stress equations in general. It is now possible to enlarge the group in order to get 
more equations, that is as many equations as the number of group parameters. Using the conformal 
group of space-time, the 4 elations give rise to 4 nonzero second order jets only which allow to 
exhibit the 4 Maxwell equations for the induction {H, D) along lines only sketched by H. Weyl in 
[20] because the needed mathematics were not available before 1970. But, as we already said, this 
is another story !. 
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REMARK 7: It becomes now clear that the ■n?{n'^ — l)/4 first order compatibihty conditions 
for the Cosserat fields [11] (the so-called torsion and curvature of E. Cartan [4]) are described by 
the second Spencer operator in the Spencer sequence while the n^(n^ — 1)/12 second order com- 
patibility conditions for the deformation tensor (the so-called Riemann curvature) are described 
by the second operator Vi in the Janet sequence. Accordingly, the torsion+ curvature of Cartan 
is not at all the generalization of the curvature of Riemann, contrary to what is still claimed in 
mathematical physics today. 

CONCLUSION: 

The original approach of V. Arnold to hydrodynamics was based on specific analytic results 
only valid for the Lie pseudogroup of volume preserving transformations. We have extended this 
approach to an arbitrary Lie pseudogroup, meanwhile showing that the previous results are in 
fact purely formal results based on a new approach to duality theory in constrained variational 
calculus. 

However, we have also found that these techniques, where the "time" variable is considered as a 
simple parameter, were not "natural" in some sense and could be superseded by the construction 
of the nonlinear Spencer sequence for an arbitrary Lie groupoid. This result, explaining the gaug- 
ing concept in the jet framework, also achieves the modern interpretation of the works done at 
the beginning of the previous century by the brothers E. and F. Cosserat on the foundation of 
continuum mechanics and by H. Weyl on the foundation of electromagnetism. 
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